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Abstract 

We study systems of few two-component fermions interacting via short-range interactions within 
a harmonic-oscihator trap. The dominant interactions, which are two-body, are organized according 
to the number of derivatives and defined in a two-body truncated model space made from a 
bound-state basis. Leading-order (LO) interactions are solved for exactly using the formalism of 
the No-Core Shell Model, whereas corrections are treated as many-body perturbations. We show 
explicitly that next-to-LO and next-to-next-to-LO interactions improve convergence as the model 
space increases. We present results at unitarity for three- and four-fermion systems, which show 
excellent agreement with the exact solution (for the three-body problem) and results obtained by 
others methods (in the four-body case). We also present results for finite scattering lengths and 
non-zero range of the interaction, including (at positive scattering length) observation of a change 
in the structure of the three-body ground state and extraction of the atom-dimer scattering length. 

PACS numbers: 03.75.Ss, 34.20.Cf, 21.60.Cs 
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I. INTRODUCTION 



The details of an interparticle potential of finite range R are largely irrelevant in few-body 
systems when the two-body scattering length 02 ^ R- This feature, which is referred to 
as universality, and its associated simplicity make such systems fertile testing ground for 
methods designed to tackle larger ones. While the measured two-nucleon scattering length 
happens to be large compared to the inverse of the pion mass, use of magnetic fields to 
create Feshbach resonances in cold, trapped atomic systems has opened up the possibility 
of dialing two-atom scattering lengths to values much larger than the typical range of the 
van der Waals potential. A recent ground-breaking achievement [H is the ability to further 
confine just a few atoms in nearly isolated sites of optical lattices formed by laser beams. At 
low temperatures, the lattice sites may be considered as harmonic oscillator (HO) traps. As 
long as the HO length b is large, b ^ R, a trapped system still displays universal properties, 
although for 6 ^ 02 these are rather different from those in a large trap or free space. 

At sufficiently low energies, an effective field theory (EFT) has been formulated which 
uses the separation between the scales 02 and R to build an expansion in powers of R/a2 
[2|-0]. It replaces the underlying interparticle potential by a series of contact interactions 
with increasing numbers of derivatives of delta functions, akin to the multipole expansion 
in classical electrodynamics. Except for isospin, which does play a role in the relative 
relevance of few-body forces, the version of this EFT used in nuclear physics [sl is formally 
indistinguishable from the one describing atomic systems [6|. (Nevertheless, the underlying 
theories for the two cases are very different.) 

As a consequence, atomic systems characterized by large scattering lengths can be studied 
with techniques developed in nuclear physics and, conversely, can provide tests for few- and 
many-body methods that can be further applied, with little or no change, to the description 
of nuclear systems at low energies. One such method is the No-Core Shell Model (NCSM), in 
which the solution to the non-relativistic Schrodinger equation for many-nucleon systems is 
obtained by numerical diagonalization using a discrete single-particle basis, typically a HO 
basis 0]. This method is characterized by its flexibility, being able to reach medium-mass 
nuclei with no limitation to closed or nearly-closed shells and, at the same time, handle local 
and non-local interactions on the same footing. The cost of diagonalization is controlled by 
restriction to a "model space" with a maximum number N"^""^ of accessible shells above the 
minimum conflguration. The NCSM is well suited to handle EFT interactions [8|, which 
are in general non-local (in the sense of involving momenta) and always deflned up to a 
maximum momentum A, the ultraviolet (UV) cutoff. 

In Ref. jsf, we formulated an EFT for two particles that support 5- wave interactions with 
02 ^ -R in a HO trap of length b ^ R. The EFT interactions are deflned within a model 
space with cutoff A = + 'i/b. We considered explicitly the flrst three orders in the 

expansion, up to which point only the S wave is modifled by the short-range potential. In 
leading order (LO), there exists only one, non-derivative contact interaction, which captures 



the physics of the scattering length and thus reproduces the results [10[ obtained with a 
pseudopotential Subleading orders involve derivatives of contact interactions treated 

in perturbation theory. The next-to-leading-order (NLO) and next-to-next-to-leading-order 
(NNLO) corrections account for the physics of the effective range r2 and lead to generalized 
pseudopotential results [l2-14|. Higher orders can be calculated similarly. We showed how 
convergence of the theory as function of N^""^ can be improved systematically, and is in fact 
better than for an exact diagonalization of subleading orders. 
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In this paper we discuss systems of three and four two-state fermions in a trap. For 
these fermions, which we can think of as having spin s = 1/2, the approach of Ref. 
apphes and we again consider exphcitly the EFT to NNLO. Thanks to the Pauh principle, 
three- and higher-body interactions appear only at higher orders [l5|, so the properties of 
few-body systems can be predicted. At LO, we solve the few-body problem using the NCSM 
formalism, as done before [3]. Beyond LO, we employ many-body perturbation theory. 

A critical aspect is the use of different values of the model-space cutoff for systems with 
different number of particles, which leads to improved convergence. This is because the 
spectator particles in a many-body system can carry some of the excitation energy, leaving 
less available to an interacting pair. If the two-body and many-body spaces are cut off at the 
same number of quanta, some states of the two-body system will simply be omitted from the 
many-body space without their effects being taken into account by renormalization. One 
way to avoid this is to use a cutoff for the space where the many-body system is solved that 
is higher than that in the two-body subsystem where the two-body interaction was defined. 
A related approach, in which the two-body cutoff is taken to depend on the state of the 



spectators, was originally considered within the NCSM formalism in Ref. [17[, but it was 
later abandoned in favor of the simpler approach with equal cutoffs. 

We calculate the lowest levels of the three- and four-fermion systems allowing for finite 
values of both 02 and r2. At unitarity, our three-fermion results converge to the semi- 



analytical values of Ref. [18[. Our improved convergence allows to better pinpoint the 
scattering length at which the ground-state parity changes, a phenomenon first identified in 
Ref. [3l and subsequently confirmed [l9|. We also investigate atom-dimer scattering and 
compare with Ref. M. In the case of four fermions, we find a ground-state energy similar 



to values in the literature |21|-|23l]. and give an example of an excited level. Systems with 
larger numbers of fermions can be dealt with at the expense of more computer time. 

Our approach has similarities to the one followed in Ref. [2l|, where an effective short- 
range interaction is fitted to several levels of the pseudopotential at unitarity, but diago- 
nalized exactly. To the extent that they are model-independent, results from finite-range 
potentials [l2[ should be equivalent to ours. 

The paper is organized as follows. In Sec. [TTl we first briefly recall how interactions 
are generated using EFT and then outline their solution in the NCSM formalism. In the 
following section. Sec. IIIH we show results for three and four fermions. We conclude and 
discuss future applications in Sec. [IVl 



II. FORMALISM 

We consider a non-relativistic system of A fermions of spin s = 1/2 and mass m in a 
HO trap of frequency u. The HO potential can be decomposed into two pieces, one acting 
on the center of mass (CM) of the A particles and one on their relative coordinates. We 
denote by r*, (pi) the position (momentum) of particle i with respect to the origin of the HO 
potential. The Hamiltonian describing the relative motion of the particles is given by 

H = Ho + Y.V,, + ..., (1) 

i<j 



with 



««-^E^p.-p,r^^m-r,r, (2) 

t<3 Kj 
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Vij being the two-fermion interaction, and ". . ." denoting three- and more-body interactions. 
In the following, the zero-point of the energy scale is such that the CM motion in the trap 
is omitted. 

The HO introduces an energy scale uj, or equivalently a length scale, the HO length 




We are interested in systems where b ^ R, the range of the force between the particles. We 
discuss in this section the construction and solution of Vij for few-particle systems. 



A. Construction of the interaction using EFT 



The interactions are constructed using EFT for the trapped few-fermion system. In free 
space the two-body interaction can be characterized by its scattering phase shifts. For small 
enough values of the on-shell momentum, k <C 1/-R, the phase shifts can be described by 
the Effective Range Expansion (ERE) |2J]. Potentials that generate the same values for 
the ERE parameters cannot be distinguished at low energies: they all lead to the same 
wavefunction for distances beyond the range of the force, r > R. This universality can be 
made manifest by taking, instead of a specific finite-range potential, an interaction expanded 
ClS db Taylor series in momentum space. 

For two-state fermions, we can use the results of Ref. [9|, as long as we interpret the 
J2i<j in Eq. ([1]) as a sum over pairs of particles in different states. In this case, an 
S'-wave interaction is allowed and dominates. In the CM of the two fermions, and expressed 
in terms of relative coordinates, the interactions considered in this paper are 



V{f', f) = Co6{f')6{r) - C2 { [V '^6{f')] 6{r) + 6{f') [v^6{f)] } 

+C4 { [V'%f')] 6{r) + 6{f') [V^(5(r)] + 2 [V'^6{f')] [V^6{r)] } 



(4) 



where Cq, C2, and C4 are parameters and ". . . " denote terms of higher orders. Since 
these interactions are singular, a regularization procedure is introduced in form of a UV 
cutoff A. This separates the short- distance physics, which is not included explicitly in the 
dynamics of this problem at low energy, and the long-distance physics, which is. In order for 
observables to be renormalization-group invariant, i.e. independent of the arbitrary cutoff, 
the parameters Cj must depend on A. 

The HO provides a natural basis on which to expand wavef unctions: its eigenfunctions 
(pni, which can be labeled by the quantum numbers n (radial) and / (orbital), have energies 
(A/'a + 3/2)0;, with 

N2 = 2n + 1. (5) 

The HO basis in turn offers a natural cutoff in the form of a maximum number of shells 
included in the basis: there exists a maximum relative momentum 



3/2, 



(6) 



with iVs""""^ 



the number of quanta of the highest energy state in the basis. The wavefunction 
of the interacting system is a superposition of HO eigenfunctions within the model space, 
determined by a solution of the Schrodinger equation with the potential (jl]). 
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In the case we are interested in, 02 S> R, S-wave interactions are enhanced by powers 
of a2 over their natural size R 0]. The LO interaction is the Co term in the potential (|4]), 
which represents the physics of the scattering length 02. The other terms are corrections, in 
particular C2 at NLO and C4 at NNLO, both of which contain the effective range r2. Higher 
5*- wave ERE terms and interactions in higher partial waves appear only beyond NNLO. To 
NNLO, then, the energies of the interacting system in waves with / > 1 are not changed 
relative to the HO energies. 

For / = 0, on the other hand, the two-body energies E2-n are determined from the solution 
of the Schrodinger equation with a superposition of HO wavefunctions up to a maximum 
jymax _ 2nmax, where Umax is the largest radial quantum number carried by states in the 
two-body basis. The solution at LO is obtained by an exact diagonalization considering 
only the Cq term in the potential (jlj). The parameter Cq{N^°'^) is fixed by imposing that 
one known energy, which we take as the ground-state energy, be reproduced at any N^"-^ . 
At NLO, the C2 term is included in first-order perturbation theory and C2{N^°-^) is fitted 
to a second known energy, which we take as the first-excited-state energy. At NNLO, C2 
is treated in second-order perturbation theory while the C4 term is included in first-order 
perturbation theory, C4^{N^°-^) being fitted to a third level, which we take as the second- 
excited-state energy. 

The energies used as input are the exact two-body energies in the limit N^"'^ — >■ c>o. 
These energies can be found by solving a simple transcendental equation involving the ERE 
parameters 



r(3/4 - E2;n/2u) ^ _^ _ r2 £^ 

r(l/4-E2;„/2u;) 2a2 2h u "' ^ ' 

In LO, we include only the 02 term; at NLO and NNLO, we include also r2. The levels not 
used as input depend on N^"'^, and for A>l/R they are predictions of the method. It can 
be shown that including more corrections to the potential improves the convergence to the 
exact solutions as the size of the model space increases. More details on the construction of 
the interactions can be found in Ref. j^. 

In this way, the two-body interaction is determined at each iV^"^ from the ERE [i.e. scat- 
tering) parameters, and can be used as input into the calculation oi A > 3 systems. In 
general, for the latter we also have to include few-body interactions. For two-state fermions, 
however, the Pauli principle requires three- and more-body contact interactions to involve 
derivatives, so that the extra fermions are placed in states with different orbital quantum 
numbers. The extra bodies and derivatives come with extra factors of R, and the corre- 
sponding interactions are suppressed [l5|. To NNLO, no few-body interaction needs to be 
included here. 



B. Many-body basis and truncation 

With Vij so constructed, we can predict the energy levels of larger systems. We need 
to solve exactly for the LO potential, and then higher-order corrections are calculated in 
perturbation theory, just as in the two-particle case [9|. 

In a shell-model approach, energy eigenstates are obtained by direct diagonalization in 
a many-body basis constructed with HO wavefunctions. There are essentially two equiv- 
alent ways to construct the basis states. In one, the basis states are Slater determinants 
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constructed from single-particle HO wavefunctions in the lab frame. This leads to the wave- 
function of the CM of the system factorizing exactly from the internal coordinates, as long 
as one performs an energy truncation of the basis states such that all states up to a given 
energy N"^°-^u:, N'^°'^ being an integer, are included. While antisymmetry is easily built into 
this approach, the dimension of the basis space grows quickly. None-the-less, this method 
can be applied efficiently to systems of more than five particles. In the second approach, 
one considers states that depend only on internal (Jacobi) coordinates and the dependence 
of each of these is described by a HO wavefunction j25|. For up to five particles, this is more 
effective than a Slater-determinant basis, but the antisymmetrization becomes increasingly 
difficult as the number of particles grows. As long as the same energy truncation is applied 
to both, the Slater-determinant and Jacobi bases give identical results for the intrinsic state 
of any system. Since in the current paper we investigate only systems with three and four 
particles, we describe in some detail the internal-coordinate approach. The novel truncation 
that we employ in the few-body system becomes more transparent in Jacobi coordinates. 

We work with Jacobi coordinates defined in terms of differences between the CM positions 
of sub-clusters within the A-body system: 
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^3 



A-l 



1 



A 



A-l 



(ri + f2 H h Ta-i) - ta 



(8) 



Using these, the HO Hamiltonian (|2]) can be expressed as 



(9) 



where is the momentum canonically conjugated to C,p- 

To illustrate the construction of the basis, let us consider three particles of spin s. In this 
case, two Jacobi coordinates ^i, ^2 are necessary to describe the internal motion. The basis 
states 

(10) 



have the spatial part constructed using HO wavefunctions with quantum numbers n, I and A/", 
C respectively, with the angular momentum coupled to L, while the spin part is constructed 
by coupling three spins s to total spin 5". In Eq. ( ITOj) . A stands for the operator that 
antisymmetrizes the three-particle wavefunction. If particles 1 and 2 are already described 
by an antisymmetric wavefunction, fulfilling the condition 



\s+l 



(11) 



then A ensures the correct behavior of the three-body state under exchange of particles 1 or 
2 and 3. The basis states thus constructed are eigenstates of the unperturbed Hamiltonian 
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Hq and the energy of each three-body state can be written as (A^3 + 3)w where the quantum 
number N3 is defined by 



N3 = 2n + l + 2Af + C. (12) 

Details of the construction of a fully antisymmetrized basis from the states (fTOj) . the calcula- 
tion of two-body matrix elements of the interaction Vij in this basis, and the generalization 



to more particles and three-body forces can be found in Ref. [25[ . 

Because the NCSM is based on a direct numerical diagonalization, the basis must be 
truncated to a computationally tractable size. For a system of A particles, the model space 
is truncated by introducing a cutoff TV™'^^ defined as the largest number of quanta in the 
eigenstates of Hq used to construct the A-hodj basis. Again using the three-body system 
as a concrete example, truncating the three-body basis at iV™'^^ means keeping only states 
with HO energies such that N3 = 2n + I + 2Af + C < N^""^. 

After the truncation of the many-body space, characterized by iV^"^', the natural question 
is what would be the consistent two-body interaction. Since the two-body interaction is 
defined in the two-body system, it is characterized by N^°'^ . In the conventional NCSM 
approach, it is customary to choose the truncation in the two-body system so that the many- 
body space is the minimal required to include completely the two-body space. For example, 
if we consider just 5- wave interactions, N^"-^ = N'^"'^ when one describes positive-parity 
states {N^°-^ has to be even), and iV^°^ = N^°-^ — 1 for negative-parity solutions (iV™"^ has 
to be odd). This was also the procedure used in Refs. [sl, [l6| . 

However, one has to consider that the renormalization of the two-body system means that 
states lying above the cutoff N^°-^ = 2nmax have been "integrated out" rather than simply 
discarded. Their effects are thus included implicitly in the effective two-body interaction. 
When these two interacting particles are embedded in a system with a larger number of 
particles, the spectators will carry energies associated with the HO levels they occupy. For 
example, of the {N^ + ?))uj total energy of one of the basis states (ITU]) . (2A/' + C + 3/2)a; 
is carried by the relative motion of the spectator. As a consequence, the maximum energy 
available to the two-body subsystem is smaller than that allowed by the A-body cutoff 
jymax some of the states removed by the truncation will not be accounted for by the 
renormalization. One way to correct for this is to use the interactions renormalized with 
a state-dependent two-body cutoff N^°-^ = N^"-^ — (2A/' + £), as first suggested within an 



NCSM approach in Ref. [17[. However, the resulting state-dependent interaction is difficult 
to handle in Jacobi coordinates for systems with more than three particles, and cannot be 
incorporated in a Slater-determinant basis. In order to account for all the two-body physics 
beyond our cutoff without the use of such an interaction, we simply decouple the cutoff 
of the many-body problem from that of the subcluster defining any interaction. Such a 
prescription has some similarity to the truncation used in Ref. [2l|. 

Each of our calculations is characterized by two cutoff parameters: A^™"^ for the two- 
body subsystem, and iV™"^ for the few-body system. To the order we work, no three-body 
forces appear and so we do not need to consider a separate cutoff for renormalization of 
a three-body subsystem, when considering larger systems. Since to this order one has to 
include only S'-wave interactions, N^"-^ is even, and N^°'^ is even (odd) for few-body states 
with even (odd) L. Our final results are obtained by first increasing 7V™"^ at fixed Nl^°-^ 
until they converge, and then increasing 7V^"^. For two-body states with > N^"-^, we 
simply set the interaction matrix elements to zero. As we increase N^"-^ (at fixed N^°'^) 
from either A^™"^ = ]\j^o-x j\jm.ax _ j\jmax _|_ observe a rapid dependence on N^°-^ 
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until it is somewhat larger than N^"'^, the difference reffecting the typical number of quanta 
carried by the spectators. For low-lying many-body states, further enlarging N^"'^ makes 
little difference because we are adding only two-body states where the two-body potential 
is switched off. Having achieved results for any observable of interest that are stable with 
respect to N^"-^ for each N^"-^, we can then take the limit of those values for large N^"-^. 
Examples are given next. 



III. FEW FERMIONS: RESULTS 

In this section we present explicit results for energies of systems made of a few two-state 
fermions in a harmonic trap. Our goal is to show convergence as we increase the UV cutoff, 
iV™"^, and its systematic improvement as the order in the EFT increases. 

Nothing in our method is specifically tied to unitarity {b/a2 = and r2/b = 0) and we can 
carry out calculations for finite scattering length, as well as finite range, as long as 02 S> ^2. 
We can in fact obtain the b/a2 and r2/b dependences of any state, improving on the r2/b = 



results of Ref. [16| because of both the accelerated convergence stemming from subleading 
orders and the inclusion of the finite range of the interaction (which in the absence of r2 
arises from the truncation of the two-body space). In this first study of subleading orders we 
limit ourselves to three and four fermions, where we can more extensively confront existing 
results. Systems of more particles can be attacked with larger computer resources. 



A. A = 3 system 

We consider first a system composed of three fermions, which allows us to illustrate the 
basic ideas of the method. Before considering the more realistic cases of finite scattering 
length and small effective range, we start with the unitary case, where semi-analytical re- 
sults exist [l8| . The three- fermion spectrum offers an excellent testing ground for numerical 
approaches, including the LO considered in our previous publication [l6| and other meth- 



ods to construct effective interactions (21|. We also use three-body energies to show how 



scattering information can be extracted for trapped few-body systems. 



1. Unitary case 

The ground state of the three-fermion system at b/a2 = and r2/b = has L = 1 and 
negative parity. Figure [1] shows the convergence of the energy of this state with the size of 
the three-body model space, ]\f^°-^, for two values of the UV cutoff in the two-body system, 
(a) N^'^^ = 10 and (b) N^°-^ = 18. For fixed N^"-^, the Hamiltonian does not change, 
so that one expects a variational behavior of the ground-state energy with increasing the 
three-body model space. Indeed, as shown in Fig. [H the energy decreases until convergence 
is reached for a large enough three-body model space. The value of N^""^ for which the 
convergence is obtained depends on the particular value of N^°-^ . Thus, for N^"-^ = 10 the 
energy of the three-body ground state (in LO and corrections) does not change by more 
than 10~^ once N^"'^ > 19, while for N^"-^ = 18 convergence at this level is achieved for 

Even though for fixed N^"-^ the errors induced by the three-body cutoff are eliminated, 
the errors induced by the truncation in the two-body sector, where the interaction is defined. 
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FIG. 1: Energy in units of the HO frequency, E/uj, of the ground state = 1~ of the A = 3 
system at unitarity, as function of the three-body model-space size, N^'^^: (a) 7V|""^ = 10; (b) 
j\j-max _ -|^g^ (Black) Circles correspond to LO, (red) squares to NLO, and (green) diamonds to 
NNLO. The (black) dashed line marks the exact value 18 1. 



can be eliminated either by taking 7V™"^ to large values or by adding corrections that take 
into account physics left out by the truncation to a certain order, or by combination of the 
two. Figure [2] shows the convergence with respect to N^°'^ for the ground-state energy at 
unitarity. The LO calculation converges to the exact result [l8|, as has been shown before in 
the case N^°-^ = N^°-^ + 1 1^. However, in Ref. 16 1 the ground-state energy had a faster 



running to the exact value: since convergence in N^"'^ is from below but in N^"-^ from above, 
increasing N^""^ at fixed N^""^ pushes the ground-state energy further away from the exact 



result. For example, when iV™"^ = 22 the LO ground state is ai E /uj = 2.7413 in Ref. (16 
while the present approach gives E/u = 2.7386. On the other hand, adding corrections to 
the potential speeds up the convergence: at NLO the agreement with the exact calculation 
is achieved faster than at LO, and improves still at NNLO. For the same N^""^ = 22, the 
energy at NNLO is E/u = 2.7715, very close to the exact solution E/u = 2.7727 found in 



Ref. [18|. Thus, subleading orders provide significant improvement over LO results. 

The calculation of A = 3 excited-state energies can be carried out similarly. We show in 
Fig. E] the running with the three-body cutoff of the energy of the first excited state with 
L'^ = 1~, for the same two values of N^""^ considered before, and in Fig. |4]the convergence 
with N^""^ . The same 10~^ precision for the same two-body UV cutoffs considered before 
is achieved at somewhat larger three-body cutoffs, N^'^^ > 23 and N^'^^ > 35 respectively. 
Like for the ground state, for a fixed N^""^ the values of energies at all orders decrease until 
convergence is reached. Note the sharp decrease of the energy as N^""^ goes from N^""^ + 1 
to N^'^^ + 3, followed by small change as the three-body cutoff is further increased. This 
suggests that a small number of quanta is carried out by the spectator, so that most of 
the two-body physics can be accommodated by a relatively small three-body space. The 
importance of having two different cutoffs in the two- and many-body systems clearly appears 
in this case. Indeed, if N^"-^ is fixed at NJ^^"-^ + 1 one can see from Fig. [3] that as corrections 
to the potential are added, results get worse: for both values of N2 the energy at NLO and 
NNLO is farther away from the exact value than the value obtained at LO. As Fig. H] shows, 
once N^"-^ is decoupled from iV™"^, the corrections to the potential again improve the energy 
systematically (except at very low two-body cutoff, i.e. N^"-^ ^10). Agreement with the 
exact value ^ E/u = 4.7727 is very good: for N^""" = 22, we find at LO E/u = 4.7457, 
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FIG. 2: Energy in units of the HO frequency, E/uj, of the ground state L'^ = 1 of the ^ = 3 
system at unitarity, as function of the two-body cutoff, N^'^^ . Notation as in Fig. [TJ 
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FIG. 3: Same as Fig. [H but for the first excited state with L'^ = 1" 



slightly below the value E/u = 4.8554 in Ref. [16|; at NNLO, E/uj = 4.7721. 

The general features of convergence with N^""^ and N^""^ are also shown by states in 
other channels, although details vary. As an illustration, in Fig. |5]we show the energies for 
the first two states with = 0~^. In this case, the convergence with N^"-^ is not always 



18| E/io 



from below. For the lowest state, agreement with exact value 
good already at NLO: for N^""" = 26, E/u = 3.1652, 
NNLO, for the same N^"'^, the result is very close, but slightly worse, E/u 
the excited state, NLO is not as good, and at NNLO E/u = 5.1614, within 0.1% of the 



3.1662 is very 
a difference of less than 0.05%. At 

3.1641. For 



exact value (l8|], E/u = 5.1662. 



Overall, there is systematic improvement as 
inclusion of higher-order interactions. 



increases, which is accelerated by the 
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FIG. 4: Same as Fig. [21 but for the first excited state with = 1" 
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FIG. 5: Same as Fig. [21 but for the ground state (left panel) and first excited state (right panel) 
with L'^ = 0+. 



2. Error analysis at unitarity 

Truncation of the two-body space at N^"-^ corresponds to imposing a two-body UV cutoff 
at the momentum scale The errors induced by this are expected to run as inverse powers 
of tlie cutoff scale, with a leading term of order (see, for example, Ref. [26[). Since 
the terms in the effective potential represent short-range physics and thus must be analytic 
in (p/A)^, only odd powers of should appear. We therefore fit our cutoff-dependent 
energies with the form 



max 1 
2 / 



E.(oo) 



ai 



as 



"5 



(^j^max 3/2)1/2 (^^max ^ 3/2)3/2 (^f^max ^ 3/2)5/2 



(13) 



expressed, for convenience, in terms of dimensionless quantities. Here £'3(00) is the semi- 
analytical result of Ref. 18|| for the energy of the state. Some care is needed with this 
procedure since, with data for only a limited range of values of N^^^, the fits can become 
unstable if too many terms are included. 

The ground-state energy at LO in the range N^"-^ = 18 to 30 can be well described by 
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FIG. 6: Energy in units of the HO frequency, Ejio^ of the lowest states with Li^ = 1~ (circles) 
and L'^ = 0^ (squares) at NNLO, as a function of the ratio 6/02 between the HO length and the 
two-body scattering length, when r2/b = 0. At 6/02 — 1-34 there is an inversion of ground state. 

fits with leading coefficients ai = —0.139 and 0^3 = —0.63. At NNLO, fits to the ground 
state lead to values of ai that are consistent with zero. Setting «i to zero and refitting the 
NNLO energy over same range of N^"-^ yields = —0.30 and as = —5.8. The greatly 
improved convergence at NNLO that is obvious from Figs. 2 and 4 is due to elimination of 
the leading (A~^) term and reduction of the coefficient of the next one (A~^). 

The coefficients of the higher-order terms in these expansions are large, implying that the 
series in A~^ is rather slowly converging for the values of the cutoff used here. This problem 
appears to be worse for the excited states, where we have not been able to find stable fits 
[i.e. fits where the values of the coefficients do not change appreciably as the number of 
terms is increased) to the energies for N^^"-^ < 30. However one should remember that the 
absolute errors on the NNLO energies are extremely small for N^""^ > 20 and so this has 
no practical consequences for our results. 



3. Finite scattering length 

Away from unitarity, universality can take on a different character. Results for the ground 
and first excited state at NNLO as a function of 6/02, but keeping r2/b = 0, are shown in 
Fig. [6l For each value of 6/02, calculations were carried out with N^"-^ = 22 and the value 
jymax increased to reach convergence (at the same level of 10~^ as before) of the energy. 
The values at 6/02 = are the same as in the previous section. As pointed out in Ref. ly 



and reproduced in Ref. [19|], at positive b/a2 there is an inversion of parity for the ground 
state from = 1~, as expected from the HO levels with a weak interaction, to L'^ = O""", 
as expected from one particle moving in an S-wave around a bound state of the other two. 
Here our use of higher orders in perturbation theory allows us to pinpoint this inversion to 



b/a2 — 1.34, very close to our original estimate of 1.5 [16[. This ground-state inversion is 
universal in the sense that it holds for all potentials with negligible r2. It is one example of 
how a trap with b^a2 can have different universal properties than a large, 6 ^ 02, trap. 
Positive scattering length means that the two-body ground state is bound as the trap 
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is removed, b — )■ oo. In this limit one is often interested in scattering on the bound state. 
Scattering observables are notoriousl y d ifficult to calculate in many-body systems, despite 



significant progress in this direction [27|. The presence of the trap, however, can simplify 
the calculation of low-energy scattering observables. In the two-body case, given energy 
levels E2-n "we can use Eq. ([7]) to extract S-wave scattering parameters. In Ref. jof we 
examined the size of such parameters induced by our two-body truncation, and as the 
present manuscript was being written a study of the extent to which such an extraction is 
possible in the two-nucleon system within the NCSM scheme appeared j28|. Because similar 
connections between scattering properties and the spectrum of the trapped system exist for 
few-body systems, we can use our three-body energies to extract parameters for scattering 
of a particle on the two-body ground state. A similar procedure to determine the scattering 
length from the full three-body solution was presented in Ref. [29| . 

As noted in Ref. [l6[, in the limit b ^ a2, the lowest three-body energy approaches the LO 
dimer energy, — l/2/xa|, which corresponds to the threshold for the scattering of one particle 
on the bound state of the other two. Indeed, if one allows the dimer to form inside a wide- 
enough trap, b > a2, the low-lying three-body spectrum can be associated with the spectrum 
of two particles (one composite) in a trap. This spectrum is connected to the atom-dimer 



low-energy parameters (scattering length a^d, effective range r^d, etc. ) Isl, [Tol, [iMi 

r(3/4 - {Es-n - E2-o)/2u) b' Tad ^3;n ' ^2;0 



r(l/4 - (Es;. - E2;o)/2u;) 2aad 2b' u 



+ ... (14) 



Here, E^-n — E2-o is the energy of the three-body system above the dimer ground state, while 
b' = l/yZ/Vi^ is the HO parameter length calculated with the atom-dimer reduced mass, 
fiad = 2m/3. Note that Eq. (fT4l) is valid only for atom-dimer relative momenta smaller than 
l/a2, since 02 is of the order of the size of the atom-dimer system, thus limiting the number 
of three-body states that can be used reliably for the extraction of atom-dimer properties. 

In Fig. [7]we plot the atom-dimer scattering length aad obtained at fixed b/a2 = 3 (circles) 
and 6/02 = 4 (squares). For each two-body cutoff, we assume that the shape and higher- 
order parameters can be neglected and, using the lowest two L = states obtained in NNLO, 
eliminate the atom-dimer effective range to obtain the scattering length. As shown in Fig. 
[71 the running is slower for 6/a2 = 4 than for 6/02 = 3 because the interaction is stronger 
with respect to the HO strength. One thus would like smaller b/a2 ratios so that better 
precision is reached for smaller cutoffs; however, for too small b/a2 ratios the dimer might 
not be able to form inside the trap. In the limit of large N2^°'^ our results approach the 



continuum value of Ref. 20 



4-. Finite effective range 

As one moves away from a Feshbach resonance, effects of the interaction range should be- 
come more pronounced and universality reduced to some degree. Assuming a finite effective 
range, we can predict the changes in the properties of the three-body system. In particular, 
we can investigate the effect on the position of the crossing between the lowest L = 1 and 
L = states. 

In Fig. |8]we plot the corresponding energies as a function of b/a2, for r2/b = 0.1. As 
one can see, such an r2 changes each of the energies by less than u in the interval displayed. 
One consequence is a change in the position where the L = becomes lower in energy than 
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Petrov et. al. [19] 
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FIG. 7: Atom-dimer scattering length in units of the two-body scattering length, (ladl^.^^ as a 
function of the two-body cutoff, N2^°'^ . NNLO results for 6/00 = 3 (black circles) and b/a2 = 4 



(red squares) are compared with the value (blue line) of Ref. 2C | 




FIG. 8: Same as Fig. [Gj but for r2/b = 0.1. In this case, the crossing point between the L = 1~ 
and L = 0"'" lowest states moves to 6/02 ^ 1.75, a larger value than in the absence of effective 
range. 

the L = 1 state, which is now at 6/02 — 1.75. Indeed, the effect of the finite positive range 
is a shift of the crossing point to larger values of b/a2- 



B. A = A system 



We now consider the system made of four two-component fermions in a trap. As before, 
we fix the value of the two-body cutoff N^"-^ and increase the size of the many-body model 
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FIG. 9: Energy in units of the HO frequency, E/uj, of the ground state L'" = 0"*" of the A = A 
system at unitarity, as function of the four-body model space size, N^""^: (a) N2^'^^ = 4; (b) 
jymax ^ g_ Notation as in Fig. [2 



space, defined here as A^™"^. We show only results at unitarity, although finite scattering 
length and effective range can be entertained as well. Of course, because of the larger number 
of particles we limit ourselves to smaller spaces. 

For states with zero total angular momentum, the smallest value for the four-body cutoff 
is N^"-^ = N^°-^ . Results for the convergence of the ground-state energy at unitarity with 
respect to N^°-^ are plotted in Fig. [9] for two values of the two-body cutoff, (a) N^"-^ = 4 and 
(b) N^°-^ = 8. As for the A = 3 ground state, at each order convergence is from above. For 
the two values of N^°-^ displayed, we notice the same rather sharp decrease of the energy as 
for the A = 3, L = 1~ first-excited state shown in Fig. [31 As in the latter case, improvement 
with order is visible only after this sharp decrease. Here near convergence is achieved when 
^max reaches N^""" + 4 for Nl^"'' = 4 and N^""" + 8 for N^'''' = 8, suggesting that more 
energy is taken away from the two-body subsystems than in the three-body case. 

Figure [10] shows the convergence in for the ground and first-excited states at uni- 

tarity. For each cutoff N^°'^ , the four-body model space was increased until convergence. In 
LO, the ground-state energy for = 10 is E'/w = 3.64, to be compared with E/u = 4.01 

m Ref. m. With corrections up to NNLO E/u = 3.52, which is in good agreement with pre- 
vious calculations where the ground-state energy was found to be 3.6 ±0.1 [2I], 3.551 ±0.009 



23l |. and 3.545 ± 0.003 (21|. Like in the case of three particles, improvement is significant 



and systematic. 



IV. CONCLUSIONS AND OUTLOOK 

We have considered systems made out of a few two-component fermions in a HO trap 
of length h using interactions generated by the application of effective field theory up to 
NNLO. To this order, the interactions are purely two-body and determined by the two-body 
scattering length 02 and effective range r2- 

Calculations at LO were performed by solving the many-body Schrodinger equation via 
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FIG. 10: Same as Fig. [2l but for the ground state (left panel) and first excited state (right panel) 
with L"" = 0"*" for the A = 4 system at unitarity. Notation as in Fig. [TJ 



a direct diagonalization, similar to the NCSM approach, whereas higher-order corrections 
were treated as perturbations. We have seen that, as in the two-body case jol, convergence 
is sped up by adding corrections to the potential. We have shown the necessity of using 
different values for the two- and many-body UV cutoffs in order to allow for enough two- 
body states in the many-body environment to match the two-body physics included in the 
construction of the interaction. 

By doing so, results at unitarity for the three- and four-fermion systems agree very well 
with other solutions, either semi-analytical [l8[ or using other numerical methods jUj- 
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We also have presented results for finite values of 02 and r2. We were able to more precisely 
determine the ratio 6/02 where, with vanishing effective range, there is a parity inversion 



in the ground state of the three-body system [16|. If the range is known, its effect on the 
location of this point can be calculated. On either side of this transition the ground state is 
universal in the sense of being the same regardless of the details of the potential, as long as 
its range R <^ 02- However, the existence of the transition shows that the trap can lead to 
qualitatively different behavior compared to the free system. 

Above the transition point, the ground state is the one expected from one atom moving 
in a S'-wave around a dimer. We were able to use our calculated three-body energy levels to 
obtain an estimate for the atom-dimer scattering length comparable to the value found in 
Ref. |20|. 

This work can be extended in various directions. First, as more data on trapped few- 
fermion systems [l| appear, one could determine r2 for specific atoms and predict its effects 
on the few-body dynamics. Second, one can apply the same method to other systems, such 
as bosons or fermions with more components. In these cases a three-body force appears 
already at LO jsf, whose parameter needs to be determined from the three-body system 
itself. The corresponding UV limit-cycle behavior jsf is expected to survive the presence 
of the trap. A system of particular interest is the atomic nucleus, where r2 is known and 
the LO three-nucleon force can be determined either from the triton binding energy jsj or 
from the neutron-deuteron scattering length (through the lowest energy levels of the trapped 
system of two neutrons and one proton). We can now calculate such systems systematically 
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to high orders. 
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